Nonlinear cascades in two-dimensional turbulent magnetoconvection 
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The dynamics of spectral transport in two-dimensional turbulent convection of electrically con- 
ducting fluids is studied by means of direct numerical simulations (DNS) in the frame of the magne- 
tohydrodynamic (MHD) Boussinesq approximation. The system performs quasi-oscillations between 
two different regimes of small-scale turbulence: one dominated by nonlinear MHD interactions, the 
other governed by buoyancy forces. The self-excited change of turbulent states is reported here for 
the first time. The process is controlled by the ideal invariant cross- helicity, H'^ = dSv ■ b. The 
observations are explained by the interplay of convective driving with the nonlinear spectral transfer 
of total MHD energy and cross-helicity. 
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Turbulent convection of an electrically conducting fluid 
is of major importance for the dynamics of stellar con- 
vection zones and the evolution of magnetic fields in 
these regions. It is thus necessary to better under- 
stand its nonlinear dynamics which can be strongly in- 
fluenced by advected temperature fluctuations and self- 
organization processes associated with the magnetic field. 
To this end the two-dimensional magnetohydrodynamic 
(MHD) Boussinesq approximation is applied as a simpli- 
fied model that asymptotically describes the behavior of 
plasmas under the influence of a strong mean magnetic 
fleld directed perpendicular to the direction of gravity 
— a configuration for example realized in the filamen- 
tary sunspot penumbra [l| where the magnetic field is 
nearly horizontal and in laboratory plasma experiments. 
Additionally, convection represents a natural and more 
realistic way of sustaining turbulence when studying its 
inherent properties by numerical simulations in contrast 
to the somewhat artificial forcing mechanisms that are of- 
ten applied in such studies. Buoyancy can turn the tem- 
perature field into an active scalar and, as a consequence, 
substantially modify the nonlinear spectral transport of 
energy and other ideal invariants by the respective turbu- 
lent cascade. Some progress has been made with regard 
to the understanding of homogeneous hydrodynamic tur- 
bulent convection, e.g. BSIl&SSlSBIidI- However, 
studies of convection in magnetofiuids (see, e.g., for a 
recent review) which focus on the small-scale properties 
of homogeneously turbulent states remain scarce. 

This Letter deals with an investigation of the inertial- 
range dynamics of convectively driven two-dimensional 
MHD turbulence by direct numerical simulation (DNS). 
The system displays quasi-oscillatory changes between a 
state where turbulence is dominated by buoyant forces 
and a state governed by nonlinear MHD interactions. 
The constantly changing cross-helicity, = Jg dSv ■ b, 
is the cause for this behavior since it regulates the im- 
portance of MHD nonlinearities compared to buoyancy 
effects. The nonlinear interaction of the turbulent cas- 
cades of energy and cross-helicity which leads to the self- 



excited alternation of turbulent regimes is reported here 
for the first time and yields further insight into the yet 
not fully understood dynamics of turbulent fiows. 

The system is described by the Boussinesq MHD equa- 
tions, see for example In two dimensions with x 
denoting the horizontal and z the vertical direction they 
read 
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where 9 denotes temperature fluctuations about a mean 
gradient, v is the velocity of the flow, b is the magnetic 
fleld, Lu — By ■ {V X v) stands for the vorticity and ip for 
the scalar magnetic potential, b = VV' x By. The symbol 
Gj, denotes the unit normal vector of the two-dimensional 
plane. The current density is given by j = —Aip. The 
dissipation coefficients rj, k are dimensionless kine- 
matic viscosity, resistivity and thermal diffusivity, respec- 
tively. The equations are given in non-dimensional form 
using a normalization to the time characteristic of large- 
scale buoyant motions = {ag\VTo\)^^^^ and the tem- 
perature gradient scale Lq = T^/\\7To\. Here, a is the 
coefficient of thermal expansion, g is the gravitational 
acceleration acting in the negative z-direction, and To{x) 
is the mean linear temperature profile. The relevant def- 
inition for T* is the root-mean-square (RMS) value of 
temperature fiuctuations, see e.g. l^l- In contrast to 



the usual setup, VTq is in the horizontal direction in 
order to eliminate the elevator instability It other- 
wise appears due to periodic boundary conditions in the 
z-direction and leads to the formation of coherent verti- 
cal jets that significantly degrade the quality of turbu- 
lence statistics. A horizontal VTq leads to a large-scale 
x-dependent buoyancy force that can not be balanced 
by the pressure gradient. Since this work concentrates 



2 



on the small-scale dynamics of turbulence the generation 
of large-scale vorticity by this effect is neglected in Eq. 
dll). The mean temperature gradient drives the turbu- 
lent flow (right-hand-side (RHS) of Eq. first term) 
by temperature fluctuations which couple with velocity 
fluctuations through the buoyancy force (RHS of Eq. ([T]), 
first term). The absence of magnetic dynamo action in 



2D MHD [13|, [3| would necessitate explicit driving of 
magnetic field fluctuations. However as the flow is ho- 
mogeneously turbulent and thus does not concentrate the 
magnetic field on the boundaries of convection cells, the 
decay of magnetic energy is a slow process on the resis- 
tive time scale = ^o/*? ^ ^b- Thus, this quantity 
is quasi-stationary during the simulation which extends 
over 15 large-scale buoyancy times ^b- 

The set of equations ([1])-(111) is solved on a 27r- 
biperiodic square using a standard pseudospectral 
method with dealiasing according to the 2/3 rule [ist . 
A simulation of the above-described setup is conducted 
with a resolution of 2048^ collocation points. The ini- 
tial state consists of randomly generated fields. The 
parameters of the run are set to ly — rj — 7 x 10^"' 
and K = 1.3 X 10^^ which corresponds to a Prandtl 
number Pr = v/k w 5.4 and a magnetic Prandtl num- 
ber Prm = ly/r] = 1. The nominal Rayleigh number 
is Ra « 2 X 10^. We note, however, that the value 
of Ra characteristic of Rayleigh-Benard systems with 
impermeable vertical boundaries is only of limited sig- 
nificance for the present periodic setup. All param- 
eters are chosen such as to obtain comparable wave- 
number intervals for the inertial ranges of the turbu- 
lent fields. The inherent anisotropy of the flow mainly 
affects the largest-scale fluctuations with wavenumbers 
fc < 4, while at smaller scales anisotropy becomes negli- 
gible. The total energy is given hy E — E"^'^'^ + E^ with 
^mhd = + E"^ = 1/2 JsdSiv^ + 9) and the tem- 
perature "energy" E^ — 1/2 JgdS9^. Analogously, the 
symbol e = emhd + ee ^ fg dS [{^luj^ + r]f) + k(V6')2] 
denotes the sum of dissipation rates of total magnetohy- 
drodynamic energy and temperature energy. The turbu- 
lent state which develops quickly after the onset of con- 
vective instability is initially characterized by E^ « 5.0, 
i?™ Ki 6.0, and E^ w 1.9. The values increase slowly by 
about 40% over the simulation period while staying in 
roughly constant ratios to each other. 

In the observed energy spectra the signatures of 
two distinct turbulent regimes appear which quasi- 
periodically alternate. This happens in a continuous 
transition where either one or the other signature be- 
comes more dominant. The time characteristic of the 
alternations and the duration of the clearest manifesta- 
tion of the different states is of the order of th- The sim- 
ulation period consists of two groups of time-intervals 
labeled "BO" and "IK". Angle-integrated energy spec- 
tra, Ek = J d^k'S{\k'\ - k)E{'k'), obtained by time- 
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FIG. 1: Angle-integrated spectra of normalized MHD energy, 
£,mhd ^ £;^hd/(-^^g^j^^)i/2^ compensated by fc^/^ (solid line), 

and temperature "energy", E^, compensated by k^^^ (dash- 
dotted line), time-averaged over the IK-intervals. The dashed 
line indicates fc~^/^-scaling with regard to -B™'^'*, the dotted 
horizontals mark scaling law prefactors. 



averaging over the group of IK intervals are shown in 
Fig. [H The spectrum of total MHD energy denoted by 
the solid line exhibits Iroshnikov-Kraichnan (IK) behav- 



ior, i;^hd 



IgI . 17 1, in the inertial range of scales, 
5 < fc < 20. The spectrum is normahzed by (emhd^o)^^^ 
with bo being the RMS magnetic field and it displays a 
constant prefactor Cjk ~ 1.8 that agrees well with pre- 
vious work on 2D MHD turbulence, see, e.g., [l^. The 
spectra of E'' and i?'" (not shown) exhibit comparable 
amplitudes and scaling consistent with the IK picture. 
Although the validity of this 2D MHD turbulence phe- 
nomenology is still a matter of controversy due to its 
neglect of anisotropy caused by the magnetic field, DNS 



18l . [19| appear to 



of two-dimensional MHD turbulence 
agree well with it. 

The temperature energy spectrum E^ displays a some- 
what shorter scaling range, 5 ^ fc 12, with an exponent 
of —5/3. The Obukhov-Corrsin scaling suggests that 
temperature fluctuations are passively advected by the 
velocity field. Note that this does not necessarily imply 
the same scaling for kinetic energy which in fact scales 
closer to k~^^'^. The spectrum of the mean square mag- 
netic potential (not shown) scales self-similarly over the 
same spectral interval as E^^'^ with IV'feP /c^''/^, con- 
sistent with magnetic field fluctuations under the influ- 
ence of large-scale driving, ji/'fep ~ E'^^^^k'^ Thus, 
during the IK phases the effects of buoyancy on the flow 
are negligible compared to the MHD nonlinearities which 
clearly dominate nonlinear dynamics. 

Angle-integrated energy spectra obtained by time- 
averaging over the "BO" intervals are depicted in Fig. [2l 
Although E^ (solid line) is not an ideal invariant of 
MHD the spectrum scales with the exponent —11/5 for 
4 ^ fc ^ 40 while for the temperature energy (dash- 
dotted line) approximately E^ ~ k~'^ is obtained. No 
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FIG. 2: Angle-integrated spectra of kinetic energy, com- 
pensated by k^^^^ (solid line), and temperature energy, E^., 
compensated by (dashed-dotted line). Both spectra are 
time-averaged over the BO-intervals of the simulation. Same 
meaning of dotted lines as in Fig. [1] 



clear scaling behavior is exhibited by i?™^*^ and E™. In 
the Bolffiano-Obukhov (BO) picture of convective tur- 
bulence 21, 23| the temperature field can drive velocity 
fluctuations. As opposed to hydrodynamic turbulence, 
spectral transfer is not characterized by the turnover time 
at scale £, ^nl ~ f-/vi, but by the scale-dependent buoy- 
ancy time, tl - {t/Bif/'^. Consequently, ~ k-'^/^ 
and E^ ~ fc~^^/^. The observed scaling of E'^ is in 
good agreement with the respective BO-law, in contrast 
to the behaviour of E^. Hydrodynamic test runs with 
similar parameters exhibit BO-behaviour for both fields. 
The spectra shown in Fig. [2] suggest that the investi- 
gated system operates in a modified buoyancy-dominated 
Bolgiano-Obukhov-like regime. The buoyancy force dom- 
inates nonlinear MHD interactions in the inertial range, 
i.e. the temperature is an active scalar significantly in- 
fluencing the energy transfer. 

The time intervals of both turbulent regimes are corre- 
lated with the evolution of total cross-helicity H'^ as de- 
picted in Fig. [21 It is important to note that other phys- 
ical quantities, e.g. the energies or the turbulent heat 
flux, fluctuate on different time scales and are not corre- 
lated with the quasi-oscillations. A simulation run with 
lower resolution 1024'^ confirms the stability of the phe- 
nomenon for at least 100^;,. The dissipation coefficients 
of this test run are set io v — 1.5 x 10"'^, rj = 7.5 x 10"**, 
AC = 4 X 10""*. The impact of cross-helicity in a strongly 
aligned turbulent system like in this simulation with 
p = (|v • h\)/{^E^E'^Y/^ e [0.8,0.95] on the nonlinear 
dynamics of MHD turbulence is known to lead to signif- 



icant modifications of spectral energy transfer [23|, |24 1 . 
The limiting case p — 1 corresponds to an Alfvenic state 
where the nonlinear dynamics is completely switched off. 
The effect of weakening of nonlinear MHD interactions 
in high cross-helicity states of turbulence is the key in- 
gredient in following considerations. 
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FIG. 3: Time-evolution of H (lower solid line), H (upper 
solid line), the associated flux injected by convective driving 
efy^j (dash-dotted line) and the dissipation rate (dashed 
line). The intervals of IK- and BO-states are indicated. 



In contrast to plain MHD turbulence, cross-helicity in 
magnetoconvective turbulence is not an ideal invariant 
as can be seen from 77^ = Oh^dS — + rj) judS = 



-inj 



^. The dissipative term has a larger effect 
at small scales due to the appearance of spatial deriva- 
tives whereas the term e?j injects cross-helicity predom- 
inantly at large scales. The interplay between all three 
terms of the cross-helicity balance in the performed sim- 
ulation is shown in Fig. [3) The time intervals of IK tur- 
bulence can be identified with negative (upper solid 
line) whereas the time intervals of BO-like turbulence 
are correlated with positive H'^ . The term e^jj (dashed 
line) is always negative permanently injecting negative 
cross-helicity into the system. Fig. |3] shows that if the 
small-scale cross-helicity dissipation (dot-dashed line) 
becomes larger than the absolute amount of the injected 
cross-helicity, |, the system switches from the IK state 
to the BO state and vice versa. 

In the following, an explanation for the quasi- 
oscillations is proposed. It is assumed that at the begin- 
ning the system operates in the IK regime. This choice is 
not essential, but provides a convenient starting point for 
the ensuing considerations. The IK regime implies that 
inertial-range energies and energy fluxes associated with 
Elsasser variables = v ± b are approximately equal 



I.e. E+ « E- with S± = 1/4 SgdS{z 
the corresponding direct nonlinear spectral fluxes « 
TjT. Consequently, total cross-helicity, H'^ = E^ — E~ , 
and the associated nonlinear flux in the inertial range, 
T^*^ = - « £+-£-, are minimal. Nonlin- 
ear MHD interactions are dominant in the inertial-range 
(see. Fig. |4]). However, the convective driving contin- 
uously generates negative cross-helicity at the rate 
predominantly at largest-scales. The apparently delayed 
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FIG. 4: Ratio of nonlinear spectral energy fluxes of MHD 
energy T^^'^ = J^*'""" dk'{v • [-v • Vv + b • Vb - VP] + b ■ 
[V X (v X b)]}fc/ (dashed line) and cross-helicity T^^^ = - 
T~ = 1/2 J^""""" dk'{-z+ ■ [z- • Vz+ + VP] + z" • [z+ • Vz" + 
VP]}fc' (solid line) time-averaged over IK- and BO-intervals 
with the total pressure P — p+b^ /2 and {•}k denoting Fourier 

ttC 

transformation. Insets display the spectral fluxes with Tf. > 
for k < 300 and negative beyond. 



adaption of the initially depleted direct nonlinear trans- 
fer of cross-helicity leads to an accumulation of H'~^ at 
large scales. The resulting growth of the breaks the 
balance of E'^ and E~ and, consequently, weakens non- 
linear MHD interactions together with the resulting spec- 
tral flux of E"^^'^ as shown in Fig. ID Thus, the process 
gives rise to inertial-range dynamics dominated by buoy- 
ancy forces. Therefore, simultaneously with the growth 
of the dynamics of the system changes toward the 
buoyancy dominated BO-like regime of turbulence. 

In the BO-like regime, due to E'^ E^ , the corre- 
sponding spectral fluxes and T^T are different. Hence, 
the spectral flux T^'^ and cross-helicity dissipation 
are larger. Thus, the BO-like regime leads to efficient 
annihilation of the accumulated cross-helicity. As H'~' de- 
creases, the system approaches Elsasser energy equipar- 
tition, E~^ « E~, and it returns to the IK regime of 2D 
magnetoconvective turbulence. Due to the continuous 
injection of cross-helicity by large-scale convection, mag- 
netic and velocity field of the flow are highly aligned. The 
same phenomenon could not be observed in topologically 
less constrained three-dimensional simulations. 

In summary, self-excited quasi-oscillations between an 
Iroshnikov-Kraichnan regime and a Bolgiano-Obukhov- 
like regime of turbulence are observed in two-dimensional 
DNS of homogeneous MHD Boussinesq turbulence. The 
highly aligned turbulent fields, a consequence of the con- 
vective large-scale driving, give rise to a quasi-periodic 
weakening of nonlinear MHD interactions in favor of 
buoyancy effects. The resulting Bolgiano-Obukhov-like 



regime allows the removal of cross-helicity by nonlinear 
spectral transfer strengthening the MHD nonlinearities 
again. The two distinct turbulent states emerge due 
to the nonlinear interplay of the cascades of energy and 
cross-helicity, a situation of general importance to MHD 
systems constrained by the presence of a strong mean 
magnetic field as found, e.g., on the Sun. 

The authors thank M. Proctor for an important remark 
regarding the simulation setup. 
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